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^ ! Abstract 

<S 

The shifted - I expansion technique (SLET) has been developed to 

solve for the eigenvalues of Schrodinger equation in three (3D) and two 
dimensions (2D). SLET simply consists of using 1// as a perturbation 
/\ ' parameter, where I = I — j3. (5 is a suitable shift, / is the angular 

j^ ■ momentum quantum number for the 3D-case, / = \m\ for the 2D-case, 

and m is the magnetic quantum number. Unlike the shifted large - 
N expansion theory (SLNT), SLET seems to be applicable to a wider 
number of problems of significant interest in physics. 



I Introduction 

A large number of problems of physics require solutions of Schrodinger, Klein- 
Gordon, and Dirac wave equations. Exact solutions of these equations exist 
only for a handful of potentials, the fact that has culminated into the de- 
velopment of many fascinating approximation techniques. However, highly 
accurate analytical solutions of these equations for potentials of arbitrary 
coupling constants are hard to find. 

In this work we shall consider the Schrodinger equation for spherically 
(three-dimensional, 3D) and cylindrically (two-dimensional, 2D) symmetric 
potentials. Among the possible solutions of this wave equation in such sym- 
metries exist, for example, the coupling constant perturbation theory, the 
variational method, the WKB approximation, the two-step procedure [1], 
the operator method [2], the 1/N expansion technique [3], the rational frac- 
tion approach [4], the direct numerical integration [5-7], the supersymmetric 
quantum-mechanics-based methods [8], the Hill determinant method [9], the 
multiple-step recursion relations procedure [10], and the shifted large-N (or 
1/N) expansion technique (SLNT) [11-13]. 

SLNT has surpassed the other approximations in its domain of applicabil- 
ity, effectively it puts no constraint on the potential or the quantum numbers 
involved. Moreover, its analytical expressions have yielded very accurate and 
fast converging results [11-16]. SLNT has also shown that even if the results 
of physical interest are in two [12] or in three dimensions [13,14], it is advan- 
tageous to work in N-dimensions and use 1/k as a perturbation parameter. 
k = N + 21 — a, N being the number of spatial dimensions of interest, / the 
angular quantum number, and a is a suitable shift which has the meaning of 
additional degree of freedom and is responsible for speeding up the conver- 
gence of the resulting energy series. However, the scope of its applicability 
to relativistic problems, in particular, to atomic, molecular, and quarkonium 
physics, has not been adequately explored. This is, most probably, because of 
the coupling nature of the Dirac equation that causes difficulties in inflating 
the dimensions of the mentioned equation [14,15]. The Schrodinger equation 
for 3D hydrogenic state in a uniform electric field is yet another example of 
the many problems that cannot be solved using SLNT. Therefore, in spite of 
its success in the Schrodinger equation for spherically [11] and cylindrically 
[12] symmetric potentials of arbitrary coupling constants, the shifted 1/N 
expansion has to be reformatted to widen its domain of applicability. 



In this paper a new technique is introduced to solve the Schrodinger 
equation. What we shall call the shifted-/ expansion technique (SLET) is a 
reformation of SLNT and can be easily extended in a straight forward manner 
to solve for the Dirac, the Klein-Gordon, and for any other wave equation that 
can be reduced to the form of the Schrodinger equation below, Eq.(l). SLET 
simply consists of using 1// as an expansion parameter, where I — I — /3. /is 
the angular momentum quantum number for the spherically symmetric (3D) 
problems, I — \m\ for cylindrically symmetric (2D) problems, and m is the 
magnetic quantum number. (3 is a suitable shift which is mainly introduced to 
avoid the trivial case / = 0. In addition, j3 is chosen so that the next leading 
contribution to the energy eigenvalue series vanishes. This choice is physically 
motivated by requiring the agreement between the shifted-/ expansions and 
the exact analytical results for the harmonic-Oscillator and the Coulomb 
potentials both in 2D and 3D cases [11-13]. We shall not adopt the classical 
limit, / — »■ oo, to workout the leading term in the energy series, as Imbo and 
co-workers [11] have done. We shall fix / and determine tq by minimizing the 
leading term of the energy eigenvalue series [12-15]. Like SLNT, SLET is also 
a pseudoperturbative technique in the sense that it proposes a perturbation 
parameter that is not directly related to the coupling constant in the potential 
of interest. 

In Sec. II we discuss the method for spherically symmetric potentials of 
arbitrary coupling constants and give the resulting analytical expressions 
in such away that allows the reader to use them without proceeding into 
their derivation. Inspired by the considerable attention that has been given 
to systems of reduced dimensions [12,17-21], we also discuss the method for 
cylindrically symmetric potentials of arbitrary coupling constants in Sec. III. 
In Sec. IV we shall show that the analytical results of SLET for the power-law 
and logarithmic potentials coincide with those of SLNT [11]. In the same 
section we shall also demonstrate that the contribution of each term in the 
energy series, Eq.(6), of SLET is the same as the corresponding one of SLNT 
[11]. In Sec.V we shall consider a 2D electron gas in the x — y plane in the 
presence of a hydrogenic potential and a magnetic field in the z-direction. 
Therein, we compare our results with those of Zhu et al [20]. Finally, we 
conclude and remark in Sec. VI. 



II SLET for the 3D Schrodinger equation 

To formulate SLET we start with the radial part of Schrodinger equation ( 
in H = 2m = 1 units) 

~-d 2 /dr 2 + 1(1 + l)/r 2 + V(r)] *(r) = E^(r). (1) 

where V(r) is an arbitrary spherically symmetric potential. If the angular 
momentum quantum number / is shifted through the relation / = I— /3, Eq.(l) 
becomes 



dr 2 



I 2 + (20 + 1)1 + [3(1 + [3) 



V(r) ) *(r) = Em(r). (2) 



It is convenient to shift the origin of coordinate by the definition [3,11] 

x = /V 2 (r-r )/r , (3) 

and to expand about x = to obtain 
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V(r ) + V'(r )^ + V"(r )'-^ + V" 
We also propose the expansion of E as [15] 
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where Q is a scale to be determined below. With the above expansions Eq.(2) 
becomes 
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where 
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On the other hand, the Schrodinger equation for a one dimensional anhar- 
monic oscillator has been discussed in detail by Imbo et al. [11], which in 
turn implies that 
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Comparing Eq.(8) with Eq.(9) and equating terms of same order in / one 
obtains 



E = Q/r 2 + V(r ), 
E l = % [2(3 + 1 + (n r + l/2)w] , 

o 

E 2 = ^[P(P+l)+ ai ], 
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E 3 = — «2, 
o 
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where a,\ and a 2 are given in the appendix. Here, r is chosen to minimize 
E [13-14,18], i. e.; 



dE /dr = and d 2 E /dr 2 > 0, 



(14) 



which gives, with / = ^[Q 



l-(3=^rlV\r )/2. 



(15) 



Where f3 is fixed by the requirement that E\ vanishes, then 

(3 = -[2 + (2n r + l)iu]/4, (16) 

with 



w = 2^ + r V"{r )/V\r ). (17) 

Eq.(15) is explicit in r and is the same equation as that given by Imbo et 
al. [11] to solve for r . 

Finally, the eigenvalues are calculated by 



E = E + \[(3(l + (3) + a 1 } + ^ : 



(18) 



where n r is the radial quantum number, and / is as much as k/2 of Imbo et 
al. [11] 



III SLET for the 2D Schrodinger equation 

The radial Schrodinger equation for cylindrically symmetric potentials can 
be reduced (in % = 2m = 1 units) to the form 

(Al 2 - 1) 



-d 2 /dp 2 + 



4p 2 



V(p) 



$(p) = E*(p). 



(19) 



where p 2 = x 2 + y 2 , I = \m\, and m is the magnetic quantum number . If / 
is shifted through the relation I = I — /3, Eq.(19) becomes 



<P + (1+ j3 -1/2)(I+ j3+ 1/2) + 
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Following the procedure described in Sec. II, Eqs. (3)-(6), one obtains 
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where 
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Also \ Up can be written as [11] 
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Comparing Eq.(23) with Eq.(22) we obtain 

e = q/pI + v( Po ), 

E x = % [2(3 + (n p + l/2)w] 
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Where po is given through the condition in Eq.(14) yielding 

l-P = y/p*V'(p )/2, 



1 

2< 



P = --(n p + l/2)w, 



w = 2^3 + p V"(p )/V'(p ), 



(24) 
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(27) 

(28) 

(29) 
(30) 



and n p is the radial quantum number, n p = 0, 1, 2, • • •. Finally, the eigenval- 



ues are 



E = E 



Pi 
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+ 
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(31) 



where £0 is given in Eq.(24). 



IV Applications to the 3D case 

For the sake of comparison with the results of the shifted large - N expansion 
technique SLNT [11], the power - law, V(r) = Ar u , and the logarithmic, 
V(r) = Aln(r/b), potentials are considered. The above potentials have, 
however, been used in heavy quarkonium spectroscopy [11,15-17]. Our results 
are given in such a way that the comparison with SLNT is made clear. 



IV. 1 Power - Law Potential V(r) = Ar v 

For the power - law potential, V(r) = Ar u , the general formalism developed 
above leads to 
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It is worth mentioning that the exact energy eigenvalues for the harmonic 
- oscillator, B 2 r 2 /4, and the Coulomb, — 2/r, potentials are obtained from the 
leading term, E , in Eq.(6) as E = B(2n r + l + 3/2) and E = -l/(n r + / + l) 2 , 
respectively, where higher - order terms vanished identically. 



IV. 2 Logarithmic Potential V(r) = A\n(r/b) 

Aln(r/b), one may simply obtain the 



For the logarithmic potential, V(r 
following results 
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For both of the potentials above, it should be pointed out that each term 
in Eq.(6) has the same contribution as that of the corresponding one in 
SLNT [11]. Mathematically speaking, (E ) S let = (E )slnt, {E 2 /1 2 )slet = 
{E 2 /k 2 ) S LNT, and {E 3 /1 3 ) S let = {E 3 /k 3 ) S LNT- 

As found previously [11], Eqs.(22)-(24) along with Eq.(6) yield remark- 
ably good results even for large values of n r . It can easily be checked that the 
rate of convergence of the three terms in Eq. (6) is approximately the same as 
00 as it is for n r — ► 0. So, if the results for small n r are accurate, one 



n 



expects roughly similar accuracy for all n r . For more details on the accuracy 
of the above predictions, Eqs.(34)-(37) for the power - law and Eqs.(40)-(43) 
for the logarithmic potentials, the reader may refer to Imbo et al [11]. 
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V Application to the 2D case 

In this section we shall consider the potential that describes a 2D hydro- 
genic donor in the presence of a magnetic field of arbitrary strength ap- 
plied perpendicular to the 2D plane. If we adopt the symmetric gauge 
A=(B/2)(— y, x, 0), we can express this potential as (in cylindrical polar 
coordinates) [12,20,21] 

\/(p) = -2/p + m 7 + 7 y/4. (44) 

Where m is the magnetic quantum number and the magnetic field strength 
is given through the parameter 7; 7 ~ B. The units of energy and length are 
the effective Rydberg and the effective Bohr radius, respectively. For more 
details on the potential in Eq.(44) and the parameters involved in it, one 
may refer to references [12,17-21]. 

Considering Eq.(31) along with Eq.(24) and (28)-(30), we have obtained 
the well known limiting values of the energies of such system at the zero and 
high-magnetic-field limits [23] as 

Edonor = -(n p + \m\ + 1)~ 2 , (45) 

and 

E L andau = l{2n p + \m\ + m + 1), (46) 

respectively. Herein, we shall report that E donor and E Landau are obtained 
by the leading term E Q of Eq.(31), where higher-order terms have vanished 
identically, i.e.; E 2 = E 3 = 0. 

For the donor state in an arbitrary magnetic field, we have numerically 
solved Eq.(28) through Eqs.(29)-(30) and Eq.(24) to find these states by 
Eq.(31). 

In Fig.l (to be supplied by authors), our results (curve of long dashes 
marked with solid circles) for the Is state show excellent agreement with 
the other results in the weak-field regime. In the strong-field regime our 
results fall in between the results of the perturbation treatment [24] (curve 
of small dashes) and the direct numerical integration [23] (solid curve) on 
which Zhu et al [20] have located their predictions (solid circles). Our results 
are, however, unique in that they tend to approach the strong- and weak- 
field perturbation theory results. The perturbation theory coupling constants 
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were appropriately defined in these regimes [24]. Likewise, we believe, should 
be the tendency of the results of any approximation technique. 

Fig. 2 (to be supplied by authors) shows the ground Is and the 2p_ states 
being weakly affected by the magnetic field since the Coulomb interaction 
dominates over the magnetic interaction for low-lying states. The higher 
excited states, on the other hand, are more weakly bound and in this case the 
parabolic quantum well, which is formed by the magnetic field, determines 
the energy spectrum. Our results (solid curves marked with solid circles) for 
the Is, 2p_, and 2p + agree with those of the series expansion calculations 
[20] (solid curves). Whereas, the results for higher excited states converge 
more rapidly to Landau levels (curve of small dashes) than those of Zhu et 
al [20], especially in the strong-field regime wherein the parabolic quantum 
well is the dominating interaction that determines the energy spectrum. 



VI Conclusions and Remarks 

In this paper, the shifted - / expansion technique (SLET) has been developed 
to solve 3D and 2D Schrodinger equations. The development is a reformation 
of the shifted large - N expansion technique (SLNT) that widens its domain 
of applicability. 

For the spherically symmetric (3D) power - law and the logarithmic po- 
tentials, both SLET and SLNT have yielded the same analytical results. 
The observations (E_ ) SLET = (E ) SLNT , {E 2 /l 2 ) SLET = (E 2 /k 2 ) S LNT, and 
(E 3 /1 3 )slet = (E 3 /k 3 )sLNT are significant for the conclusion that SLET is 
indeed a reformation of SLNT. The accuracies and the speeds of conver- 
gence of both techniques are the same, at least for the power - law and the 
logarithmic potentials. 

For cylindrically symmetric (2D) potentials, we have considered a 2D 
Coulomb field in the presence of a magnetic field of arbitrary strength, 
Eq.(44). SLET results have appeared to be fast converging in the sense that 
the dominating contribution to the energy series, Eq.(6), is that of the lead- 
ing term E^-The results were compared with those of the perturbation [24], 
the direct numerical integration [23], and the series expansion [20] methods 
according to whichever was available. Figs.l and 2 show that SLET results 
seem more appropriate than the others. Fig.l shows that they approach 
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the perturbation theory results at weak- and strong- field limits, and Fig. 2 
shows that they approach the Landau levels as one goes to higher-excited 
states where the magnetic interaction is more effective in determining the 
energy spectrum. 

In the near absence of highly accurate analytical approximation methods 
to solve the Schrodinger equation, even for simple cases such as a 3D hy- 
drogenic state in a uniform electric field [25] (to be investigated in the near 
future), SLET provides highly accurate analytical expressions. Moreover, 
the difficulties associated with the application of SLNT to Dirac equation 
[14,15] should vanish for SLET. To appear elsewhere, we have applied SLET 
to Dirac, and Klein-Gordon wave equations. Finally, we should like to point 
out that the energy states of excitons in a harmonic-quantum-dot [26] , and of 
shallow donors and heavy - hole excitons in quantum well in the presence of 
magnetic field can be correctly obtained by SLET and compared with others 
[19,23]. 

Appendix 

Although some of the following definitions can be found in some other 
references we would like to repeat them so as to make this article self - con- 
tained. The definitions of a\ and a 2 and the parameters involved are as 
follows 



ai = [(1 + 2n r )e 2 + 3(1 + 2n r + 2n 2 r )e 4 ] 

- ur 1 [e2 + 6(l + 2n r )eie 3 + (11 + 30n r + 30n^], (47) 

a 2 = (1 + 2n r )d 2 + 3(1 + 2n r + 2n 2 )d 4 

+ 5(3 + Sn r + 6n 2 r + 4njJ)d 6 

- w _1 [(l + 2n r )e\ + 12(1 + 2n r + 2n 2 )e 2 e A + 2e 1 d 1 

+ 2(21 + 59n r + 51n^ + 34n^)e^ + 6(1 + 2n r )eirf 3 
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with 



+ 30(1 + 2n r + 2n 2 .)e 1 d 5 + 6(1 + 2n r )e 3 d 1 

+ 2(11 + 30n r + 30n 2 r )e 3 d 3 + 10(13 + 40n r + 42n 2 r + 28n^)e 3 rf 5 ] 

+ ur 2 [4e?e 2 + 36(1 + 2n r )e 1 e 2 e 3 + 8(11 + 30n r + 30n 2 r )e 2 e 2 3 

+ 24(1 + n r )e\e A + 8(31 + 78n r + 78n 2 )eie 3 e 4 

+ 12(57 + 189n r + 225n 2 r + 150n^)e^e 4 ] 

- w~ 3 [8ele 3 + 108(1 + 2n r )e\e\ + 48(11 + 30n r + 30n;?)eieJj 

+ 30(31 + 109n r + 141n2 + 94n^)e|], (48) 



e ^ = ^h and ^ = ^k' ( 49 ) 



where j = 1,2,3,4, i = 1, 2, 3, 4, 5, 6. 
Ej 's and Si 's for 3D case are given as 

£i = -2(2/3+1) , £ 2 = 3(2/3 + l), (50) 

E3 -~ A + ^6Q~ ' £4 " 5+ 24Q ' (51) 

<J 1 = -2/3(1 + /3) , 5 2 = 3/3(1 + /3), (52) 

5 3 = -4(2/3+1) , 5 4 = 5(2/3+1), (53) 

55 = " 6 + ^20^ ' 4 = 7+ -T20^- (M) 
For the 2D case Ej 's and Si 's are given by 

£i = -4/3 , £ 2 = 6/3, (55) 
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e - l i Pq V '"M £ -5 1 P" V ""^ (56) 

£ *-- A+ QQ ' £4 " 5+ 24Q ' (56) 

8l = - 2 {(3 2 - 1/4) , 5 2 = 3(/3 2 - 1/4), (57) 

<J 3 = -8(3 , <5 4 = 10/5, (58) 

^ = -6 + ^ j Po) , 4 = 7+ ^ M . (59) 

5 120Q ' 720g v ; 
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